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The azimuthal anisotropy of charged particles in heavy ion collisions is an important probe of 
quark-gluon plasma evolution at early stages. In the present paper the elliptic flow pattern in 
Pb+Pb collisions at y/s — 5.5 TeV is analyzed for different hadron species in the frameworks of 
HYDJET++ Monte-Carlo model. The influence of resonance decays on particle flow is investi- 
gated. The different methods of elliptic flow reconstruction are compared under LHC conditions. 
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1. Introduction 

In non-central collisions between two nuclei the beam direction and the impact parameter 
vector define a reaction plane for each event. The observed particle yield versus azimuthal angle 
with respect to the event-by-event reaction plane gives information on the early collision dynamics 
[J|, g|. An initial nuclear overlap region has an "almond" form at non-zero impact parameter. 
If the produced matter interacts and thermalizes, pressure is built up within the almond shaped 
region leading to anisotropic pressure gradients. This pressure pushes against the outside vacuum 
and the matter expands collectively. The result is an anisotropic azimuthal angle distribution of the 
detected particles. One can expand this azimuthal angle distribution in a Fourier series. The second 
coefficient of the expansion vi is called the elliptic flow. 

It was found [Q, ||] that anisotropic flow is self-quenching phenomenon since it reduces spatial 
anisotropy as it evolves. Therefore, observed elliptic flow must originate at early stages of the 
collision when the anisotropy is still present in the system. There is no elliptic flow generated 
when the spherical symmetry is restored in the system. Thus the elliptic flow keeps information 
about hot and dense matter created in relativistic heavy ion collisions. 

As the fireball expands, its temperature and energy density drop. Finally, at the freeze-out 
stage, the system breaks up into hadrons and their resonances. The effect of resonance decays, i.e. 
final state interaction, on the resulting elliptic flow of particles is quite important. For instance, it 
can explain partly the observed deviation of pion elliptic flow from the so-called constituent quark 
scaling The exact resonance decay kinematics at very low momenta can probably be accounted 
for the reduction of V2 coefficient for pions at midrapidity in the hydrodynamical calculation 

In this work we employ the HYDJET++ model [|6|, 0] to estimate the azimuthal anisotropy 
of particles in Pb+Pb collisions at LHC energy and to study the influence of resonance decays, jet 
production and jet quenching on elliptic flow, and also to test the different methods of V2 restoration 
for different particles. 

2. Simulation of elliptic flow with HYD JET++ model 

The HYDJET++ model ^ [7|] represents a superposition of soft and hard parts. These parts are 
independent and their contribution to the total multiplicity production depends on collision energy, 
centrality and is tuned by model parameters. The hard part of the model is identical to the hard part 
of the HYDJET model [^] and has the possibility to account for jet quenching effect. The soft part 
of HYDJET++ event represents the "thermal" hadronic state where multiplicities are determined 
under assumption of thermal equilibrium [[J. Hadrons are produced on the hypersurface repre- 
sented by a parametrization of relativistic hydrodynamics with given freeze-out conditions. Feed 
down of hadronic resonances is taken into account. HYDJET++ is capable of reproducing the bulk 
properties of heavy ion collisions at RHIC, i.e. hadron spectra and ratios, radial and elliptic flow, 
femtoscopic momentum correlations, as well as high-pj- hadron spectra [|6],[7]]. 

Fig. [j] presents the impact parameter b dependence of the elliptic flow coefficient V2 of charged 
particles. The coefficient V2 is defined here as the cosine of twice the azimuthal angle of a particle 
relative to the reaction plane angle (which is known in each simulated event), and averaged over 
all charged hadrons in each event, V2 = (cos(2(p — *P#)). As expected, the elliptic flow coefficient 
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Figure 1: The impact parameter dependence of elliptic flow in HYDJET++ model for Pb+Pb collisions at 
LHC energy. 
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Figure 2: The /?7- -dependence of elliptic flow in HYDJET++ model for different hadron species produced 
in Pb+Pb collisions at s/s = 5.5 TeV with centrality 42%. 



grows with increasing impact parameter (i.e. with increasing of azimuthal anisotropy of initial 
nuclear overlap region). Further we consider elliptic flow at fixed centrality, a /o geo = 42% (b 
1.3/?). Fig. § shows the pj dependence of elliptic flow coefficient for most abundant hadrons, 
i.e. pions, kaons, protons, lambdas and sigmas. This behavior of elliptic flow can be explained in 
following way. The flow of hydro part rises monotonically up to V2 ~ 0.5 at pj ~ 6 GeV/c while 
the relative contribution of hydro part to particle multiplicity decreases with pj, so the particles 
with pj > 6 GeV/c are produced only through jets (hard part). The flow of the jet part is close to 
zero in this pj range. It results in initial rise of v^{pt) followed by the fall off at pj > 3 GeV/c. 
The jet part at pj > 7 GeV/c also presents some amount of flow (« 4%) at the LHC energies due 
to jet quenching effect. The energy loss of the high pr partons depends on the passing length of 
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Figure 3: The /^-dependence of elliptic flow in HYDJET++ model produced in Pb+Pb collisions for all 
hadrons and for direct hadrons at s/s = 5.5 TeV with centrality 42%. 



the anisotropic matter thus giving the different yield of the high-p^ partons in the in-plane (x,z) 
and out-of plane (y,z) directions. 



The pronounced feature of the RHIC experimental data fllOfl reproduced by HYDJET++ in 
Fig. H is the crossing of baryon and meson branches. In HYDJET++ at low pj the hydro part 
dominates and the flow is strictly ordered by particle masses. The lighter particles (pions, kaons) 
have the larger flow than heavier (protons, lambdas). 

The slope of the heavy particles /^-spectra is steeper then one of the light particles, as a result 
the hydro part dominates till larger pj values for heavier particles: i.e. for pions hydro dominates 
till ~ 4 GeV/c, for protons it dominates till ~ 5 GeV/c. As a result at pj > 4 GeV/c the mass 
ordering changes on the opposite: the heaviest particles have the largest flow. 

The influence of resonance decays. At RHIC energies the transition from baryon rich matter 
to meson reach matter was found. As was predicted by Hagedorn, at high energies most of the 
particles will be produced through resonance decays with shifting of the average mass to heavier 
sector. The effect of resonance decays should be accounted for when one considers the V2/n q scal- 
ing. Table |l] shows the contributions of direct and resonant production for various hadron species 
including feed down from weak decays for our Pb+Pb event sample generated with HYDJET++ at 
LHC energies. One can see that 80% of pions, 70% of protons, 61% of £ and A-hyperons and 56% 
of kaons are produced from resonance decays. Figures || and |] display difference between V2 of 
all these hadrons and V2 of only direct hadrons. A degree of the influence of resonance decays on 
elliptic flow coefficient is quite different for various hadrons. The effect is strongest for protons; 
rather moderate for £ and A-hyperons and pions; and negligible for kaons. 

Let us consider the case of pion and proton flow. The relative contribution from resonance 
decays for them is presented in Table ^[ Fig. |5| shows differences in secondary pions and protons 
from (anti)deltas (A ++ , A + , A , A~) decays. Because of the kinematics of decay, when heavy baryon 
resonance decays into secondary baryon plus pion, the most part of its momentum is carried by the 
baryon while the pion is produced with low pj. Thus, the pion elliptic flow gets an extra boost 
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Figure 4: The same as Fig. 2, but for protons (upper left), pions (upper right), kaons (bottom left), lambdas 
plus sigmas (bottom right). 



Table 1: Yields of the particles produced directly and with resonance decays at midrapidity region, c=42%. 
Weak decays of the strange particles are included. 
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Table 2: Yields of the pions and protons produced directly and from resonance decays in HYDJET++, 
c=42%. 
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Figure 5: Top: The /^-dependence of elliptic flow for pions (left) and (anti)protons (right) coming from 
A decays. Bottom: The elliptic flow /^-dependence of pions coming from p decays (left) and CO decays 
(right). 



at low pj from the flow of heavy resonances (Fig. ||, top left). On the other hand, the secondary 
baryon carries practically the same V2(pr) as mother particle (Fig. ||| top right). Note, that most of 
secondary pions are produced from p- and ft)-mesons. But momentum distribution for p -meson and 
for secondary pions are quite close, thus the V2(pr) almost coincides for them (Fig ||, bottom left). 
The influence of heavy resonance decays on pions flow at pr < 2 GeV/c is effectively compensated 
by influence of p-decays. For higher pj the contribution from A'-decay determines the observed 
excess of the pion flow over direct pion flow. 

The secondary protons come from A and A decays in approximately equal proportions. They 
predominantly possess the flow of these resonances, as can be seen from top left Fig. ||. 

One can see from Fig. |5| that hadrons produced from resonance decays carry the same ampli- 
tude of V2 as that of the mother particle but the maximum may be shifted to the soft pj region. 
Elliptic flow of pions from the p — > nn decay almost coincides with v p 2 (Fig. ||, bottom left) while 
in 3-particle decay co — ► 717171 pions are getting obviously softer pj distribution, thus their elliptic 
flow is transferred to the softer pj region compared with (pr) (Fig. ||, bottom right). 

Therefore, contributions from resonances sometimes increase and sometimes decrease initial 
elliptic flow of directly produced pions. This effect, especially pronounced at LHC energy, can 
lead to violation of the mass-hierarchy in the V2(pr) sector. 
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3. Reconstruction of elliptic flow at the LHC 

There exists a wealth of anisotropic flow measurement methods, each having its own advan- 
tages and limitations. Here we apply three wide-spread methods to calculate the v 2 coefficient. 
One of them uses the event plane angle determination, and others are cumulant and Lee- Yang zero 
methods. 

The event plane angle, *P ra , can be determined from the measured ?i-th harmonics via the 



standard method [11, 12 ] : 



X>;sin(n0i) 

tannin = - -. -, n>l, < ¥ n < 27t/n, (3.1) 

X>;COs(«<p,-J 

i 

where <p,- is the azimuthal angle of the i-th particle and w, is the weight. The sum runs over all 
particles in given event. The observed value of v 2 is calculated using the event plane (EP) method 
by the formula: 

v?s{EP} = (cos2((p-*i> 2 )), (3.2) 

where event plane angle x i' 2 is the estimate of the true reaction plane angle the mean was taken 
over all charged particles in a given event and then over all events. Usually the true elliptic flow 
coefficient is evaluated by dividing by the factor R [12], which accounts for the event plane 
resolution: 

V2{EP} ~ ~ (cosier -*,)>• (33) 

This procedure relies on the assumption that there are no non-flow correlations (e.g., correlations 
due to momentum conservation, quantum statistics, resonance decays, jet production) or that they 
are negligible, and also that the full event multiplicity is large enough. Generally such assumptions 
are not true. Studying these effects in real data is a separate and non-trivial task. 

In order to avoid the trivial autocorrelation of particles the event plane angle *¥ 2 and hence R 
are calculated in one angular distribution sample of event, and v 2 in another sample with the same 
multiplicity. The samples may be selected in two regions of pseudorapidity r\ < and T] > 0. 

The basic idea of the cumulant method is that the v 2 coefficient can be expressed in terms of 
particle azimuthal correlations [13, pQ . The procedure is to construct two-particle correlator or 
cumulant 

v 2 {2} 2 = (cos2(<p i - - <pj)) = (cos2((<p ( --¥«) - 

~ <cos2((<p ; --^))<(<p;-^))). (3.4) 

As in the first method it is necessary to exclude the autocorrelations. The approximative equation 



in the last string of Eq. (3.4) means that the non-flow correlations are small. 



The two-cumulant method can be extended to the case of many-particle correlations [14]. It 



is known as the higher order cumulant method or the method with Lee- Yang zeroes. The higher 



order cumulant methods are less sensitive to non-flow effects. The Lee- Yang zeroes method [15] 
refers to the Generation function as a complex function of variable r: 

G e (ir) = (e irQe ), (3.5) 
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Figure 6: (a) The ^-dependence of elliptic flow in HYDJET++ model for Pb+Pb at LHC energies recon- 
structed by different methods. (b) The ratio of the methods to the true value. 



where Q 6 is the flow vector and B is arbitrary angle. 



<2 9 =£cos2((p ; --e). (3.6) 

j 

The integral v 2 value is connected with the first minimum of module of the Generation function 

\G°(ir)\, 

vf {-} = (3-7) 
and the differential value is given by more complex expression: 



Mv 2 {oo} ^ {Q°e<Q e ) ) 

The natural resolution parameter that appears in these methods is % ~ v%\/M, i.e. methods 
depend both on strength of the flow v 2 and on multiplicity M. The flow increases with rise of the 
impact parameter while the multiplicity decreases. The best reliability of the methods is achieved 
at midcentral collision (Fig. |6)). The overestimation of the true elliptic flow by V2{EP} and v 2 {2} 
methods in most central and peripheral collision is up to 30%. This is due to nonflow correlations. 
We see that the Lee- Yang zeroes method is good enough to reconstruct v 2 in large centrality range 
except very peripheral collisions where the multiplicity is too small. 

The origin of nonflow correlation is mostly due to jets. It can be revealed from Fig. ^ where 
the event plane method severely overestimates the original elliptic flow at high pj. The possibility 
of v 2 reconstruction for different hadron species based on the generalization of Lee- Yang zeroes 



method [16, 17] at the LHC is under investigation. 



4. Conclusion 

The elliptic flow pattern in Pb+Pb collisions at y/s = 5.5 TeV is analyzed for different hadron 
species in the frameworks of HYDJET++ Monte-Carlo model. Resonance decays and in-medium 
jet fragmentation result in the smearing of hydro-induced mass-ordering of elliptic flow coefficients 
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p T , GeV/c p T , GeV/c 

Figure 7: The /^-dependence of elliptic flow in HYDJET++ model for Pb+Pb at LHC energies recon- 
structed by Event plane method (left) and Lee-Yang zeroes method (right). Centrality is 42%, charged 
hadrons. 

V2 for different hadron species in low- and high-p^ domains, respectively. Increase of V2 due to 
resonant production is strongest for protons, moderate for £ and A-hyperons and pions, and negli- 
gible for kaons. The total effect on a given particle specie may contain two different contributions, 
either increasing or decreasing the direct flow V2(pr), as we saw on pion example. 

The comparison between three different methods of elliptic flow reconstruction under LHC 
conditions has been performed. The event plane, the two-particle correlation and Lee- Yang zeroes 
methods show the different restoration power. The event plane and two-particle correlation methods 
work well for low-pr region where jet influence is negligible, while Lee- Yang zeroes method is 
able to remove non-flow (jet) correlations at high pj. 
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